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We discuss how to build some partially entangled states of n two-state quan- 
tum systems (qubits). The optimal partially entangled state with a high degree of 
symmetry is considered to be useful for overcoming a shot noise limit of Ramsey 
spectroscopy under some decoherence. This state is invariant under permutation of 
any two qubits and inversion between the ground state |0) and an excited state |1) 
for each qubit. We show that using selective phase shifts in certain basis vectors and 
Grover's inversion about average operations, we can construct this high symmetric 
entangled state by (polynomial in n) x 2 n / 2 successive unitary transformations that 
Q-i are applied on two or three qubits. We can apply our method to build more general 

entangled states. 

• i-H 

1 Introduction 

Recently rapid progress in quantum computation and quantum information theory have 
been made|jl|[§]. In these fields, properties of quantum mechanics, which are superpo- 
sition, interference, and entanglement, are handled skillfully. After Shor's algorithm for 
factorization and discrete logarithms and Grover's algorithm for search problems appeared 
01 00 00' man y researchers have been proposing methods for the realization of quan- 
tum computation and developing quantum algorithms. On the other hand, in the fields of 
quantum information theory, it is recognized that entangled states play important roles 
for robustness against decoherence ||. 

As an application of these results, it is considered to overcome the quantum shot noise 
limit by using entangled states of n two-level systems (qubits) for Ramsey 
spectroscopy||[T(J. (M. Kitagawa et al. gave a similar idea, though an experimental 
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scheme that they discussed was not Ramsey spectroscopy of qubitsJTT|.) When we can 
neglect decoherence of the system caused by an environment, the maximally entangled 
state serves us an improvement of a frequency measurement. In this case, the fluctuation of 
frequency is decreased by 1/y/n. (In this paper, for example, we consider (l/y / 2)(|0 ■ • • 0) + 
|1 • • -1)) one of the maximally entangled states. Entanglement for n(> 3)-qubit system 
has not been defined clearly (T^j.) However, if the decoherence is considered, the maximally 
entangled state provides the same resolution that an uncorrected system provides ]TB[ . 
S. F. Huelga et al. proposed using a partially entangled state which has a high degree 
of symmetry. This state is invariant under permutation of any two qubits and inversion 
between the ground state |0) and an excited state |1) for each qubit. If we prepare the high 
symmetric partially entangled state optimized numerically, it provides high resolution in 
comparison with the maximally entangled states and uncorrelated states. 

Carrying out an experiment of Ramsey spectroscopy with the optimal high symmetric 
partially entangled state, we have to prepare it for an initial state as soon as possible, 
before time limit of decoherence. In this paper, we study how to construct this state 
efficiently We estimate time to prepare it by the number of elementary quantum gates 
that are unitary transformations applied on two or three qubits || || [[Uj] . The number 
of gates is considered to be in proportion to the amount of time for building the state. 
We show it takes 0((n 3 log 2 n) x 2 n//2 ) steps at most to build it. (It was shown that any 
unitary transformation U (G U(2 n )) can be constructed from 0(n 3 2 2n ) elementary gates 
at most||14||.) Furthermore, our method can be applied to build more general entangled 
states. 

Before discussing how to build partially entangled states, we try to construct the 
maximally entangled state with n qubits from an initial state |0 • ■ ■ 0). To do it, we need 
two unitary transformations for elementary gates. They are (the Walsh- Hadamard 
transformation) which operates on the j-th qubit, and /\i' k \cr x ) which operates on the 
j-th and k-th qubits: 
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Because /\i (cr x ) transforms \x, y) (x,y G {0, 1}) to \x,x®y) (applying a x on k-th qubit 
according to the j-th qubit), it is sometimes called the controlled-NOT gate. Applying 
Ai \ a x) ' ' ' l\ { \' 2 \o'x)H ( ^ on 1 0) i ® • • • ® |0)„, we can obtain the maximally entangled 
state, (l/ v / 2)(|0---0) + |l---l)). 

But, building partially entangled states like 



1^4) = o |0) s + a x \l) s + a 2 |2) s 
= a (|0000) + |1111>) 

+ai(|0001) + |0010) + |0100) + |1000) 

+ |1110) + 1 1101) + |1011) + |0111)) 
+a 2 (|0011) + |0101) + |0110) + |1001) + |1010) + |1100>) 
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(this is an example of the 4-qubit high symmetric partially entangled state), where do, 
ax, a 2 are given (real) coefficients, and \k) s is an equally weighted superposition of k or 
(4 — k) excited qubits, we feel difficult. It is hard to resolve a unitary transformation that 
transforms |0000) to into local operations like or /\i' k \(T x ). This is because we 
don't know a systematic method for adjusting coefficients of basis vectors. This matter 
is a motivation of this paper. 

This paper is arranged as follows. In §[| we explicitly describe the high symmetric 
partially entangled states. We make preparations for our method of building them. In 
§|3], we introduce a unitary transformation that makes two sets of basis vectors classified 
by their coefficients be weighted equally. We derive a sufficient condition for finding 
an appropriate parameter that characterizes this transformation. In §|], we develop a 
technique which transforms the state that doesn't satisfy the sufficient condition derived 
in §[3] into a state that satisfies it. This technique is an application of Grover's amplitude 
amplification process |J 0. In §||, we show the whole procedure for building the high 
symmetric entangled states and give a sketch of implementation for it. We estimate the 
whole number of elementary gates of our method. We also show that we can use our 
procedure for building more general entangled states. In §|| we give a brief discussion. 
In §Appendix, we construct networks of quantum gates for our method concretely, and 
derive a variation of coefficients of the state under the transformation discussed in 

2 High symmetric partially entangled states 

In this section, we define high symmetric partially entangled states explicitly. We also 
make preparations for our method of building them, defining an initial state, giving some 
unitary transformations used frequently, and so on. 

The partially entangled state which has a high degree of symmetry is given by 

[n/2] 

\ipn) = E a k\k)s for n > 2, (2) 

k=0 

where \n/2\ is the maximum integer that doesn't exceed n/2[0|. {a^} are given real 
coefficients. We assume a k > for k = 0, ■ • • , [n/2\ for a while. \k) s is an equally 
weighted superposition of k or {n — k) excited qubits, as shown in l^) of Eq. ([]]). This 
state has symmetric properties, invariance under permutation of any two qubits, and 
invariance under inversion between |0) and |1) for each qubit. A main aim of this paper is 
to show a procedure for building \ip n ) efficiently. We emphasize that {a k } of \ip n ) in Eq. (|2]) 
are given and numerically optimized to realize high precision for Ramsey spectroscopy. 

We make some preparations. To build \ip n ), we prepare an n-qubit register in a 
uniform superposition of 2 n binary states, (l/v^™) X^e{o,i} n \ x ) 1}" represents a set 
of all n-bit binary strings), and apply unitary transformations on the register successively. 
(Initializing the register to |0 ■ ■ - 0) and applying H^> (1 < j < n) on each qubit, we can 
obtain the uniform superposition.) 

In our method, we use two kinds of transformations. One of them is a selective phase 
shift transformation in certain basis vectors. It is given by the 2™ x 2 n diagonal matrix 
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form, 



xy ~ 1 n f or x y ' ^ 




where subscripts x, y represent the basis vectors G {0, 1}"} and < 6 X < 2-n for 

Vx. (Although a general phase shift transformation in the form of Eq. (|3|) takes a number 
of elementary gates exponential in n at most, we use only special transformations that 
need polynomial steps. This matter is discussed in §|5] and § Appendix A.) The other is 
Grover's inversion about average operation -D0. The 2 n x 2™ matrix representation of D 
is given by 

/ -l + 2-"+\ for x = y 
xy \ 2- n+1 , ioix^y ' K > 

Because we use only unitary transformations and never measure any qubits, we can 
regard our procedure for building \i/j n ) as a succession of unitary transformations. For 
simplicity, we consider a chain of transformations reversely to be a transformation from 
\ip n ) to the uniform superposition instead of it from the uniform superposition to \i[> n ). 
Fortunately, an inverse operation of the selective phase shift on certain basis vectors is 
also the phase shift, and an inverse operation of D defined in Eq. ([|) is also D. In the 
rest of this paper, because of simplicity, we describe the procedure reversely from \ip n ) 
to the uniform superposition. (Building \%p n ) actually, we carry out the inversion of the 
procedure.) 



3 Making basis vectors be weighted equally 

At first, we show how to transform \ip 2 ) to the uniform superposition as an example. After 
that, we consider of \ijj n ) for n > 3. 

Writing \tp 2 ) as 

|V> 2 ) = ao(|00) + |11)) + ai(|01) + 1 10>) where a > 0, a x > and a 2 + a\ = 1/2, 

we apply the following transformations on it. Shifting the phase of |01) by 9 and shifting 
the phase of 1 10) by (—9), we obtain 

ao(|00) + |11)) + a ie ie \0l) + a ie ~ ie |10). 

The value of 9 is considered later. Then, we apply D on the above state. D is given as 
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and we get 

A)(|00) + 1 11>) + Ai|01) + A*|10) where A = aicos9, A x = a -zaxsin^. 
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Defining <p as e 1 ^ = A\/\Ai\, we shift the phase of |01) by (—</>) and shift the phase of 1 10) 
by (f>. We get 

A (|00) + |11)) + |A|(|01> + |10)). 

If A = \Ax\, we obtain the uniform superposition. Here, we can assume < a < 
1/2 < <2i without losing generality. From these considerations, the value of 6 is given by 
cosfl = l/(2ai). 

In case of n > 3, we take the following method. Classifying basis vectors G 
{0, 1}"} of \ip n ) by their coefficients, we obtain (L^/2J + 1) sets of them characterized by 
a,k- We consider the transformation that makes two sets of basis vectors (e.g. sets of basis 
vectors with ao and ai) be weighted equally and reduces the number of sets by one. If we 
do this operation for \n/2\ times, we obtain the uniform superposition. 

Here, we consider how to make a set of basis vectors with ai be weighted equally to a 
set of them with ao on \ip n ). A similar discussion can be applied on other sets of them. 
From now, we write \ip n ) as 

1^) = K>, ■■■,cti, ■ ■ ■ ,a 2 Q+m),- •-,a 2 n-i] for n > 2. (5) 

21 2m 

As the representation of Eq. (|^), we sometimes write a column vector by a row vector. 
In Eq. (Bl), we order the orthonormal basis vectors G {0,1}™} appropriately, and 

coefficients ao and ai are put in the left side of the row. Because \if> n ) is invariant under 
inversion between |0) and |1) for each qubit, the number of basis vectors that have a 
coefficient a^ (0 < k < \n/2\) is even. Therefore, we can give the number of ao by 
21 and the number of a\ by 2m, where I > 1, m > 1, and I + m < 2 n ~ 1 . The other 
(2 n — 21 — 2m) coefficients, {a 2 , ■ • • , a\_ n / 2 \}, are gathered in the right side of the row and 
they are relabeled {aj\2(l + m) < j < 2 n — 1}. Reordering basis vectors never changes 
matrix forms of R defined in Eq. (0) and D defined in Eq. (f|), except for permutation of 
diagonal elements of R. 

We carry out the following transformations. Firstly, we shift phases of m basis vectors 
with coefficients ai by 9 and shift phases of the other m basis vectors with coefficients a\ 
by {—&). How to choose the value of 6 is discussed later. We obtain 

Re\%) = [a , • • • , e l V, • • • , e~' e ai, ■ ■ ■ , a 2 (t+m), • • • , 02»-i], (6) 
21 m m 

where < 9 < 2n (Rg is given by 2™ x 2 n diagonal matrix whose diagonal elements are 

{l,---,e*V--,e-*V--,l,---,l} )• 
Then we apply D on R g \ty), 

DR e \^) = [A , A x , ■ ■ ■ , AI, ■ ■ ■ , A 2[l+m) , • ■ ■ , A 2 n_ x \, (7) 
21 m m 

where 

2 n ~ 1 A = {2l-2 n ~ 1 )a + 2ma 1 cos6 + C, 

2 n ~ 1 A 1 = 2la + (m - 2 n ~ 1 )a 1 e ie + ma x e- ie + C, (8) 
2 n ~ l Aj = 2la + 2ma 1 cos6-2 n ~ 1 a j + C, 
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for j = 2(1 + m), ■ • ■ , 2 n — 1, and C = Y^=2a+ m ) a j- We n °tice that Ai = Aj, if = ctj for 
2(1 + m) < V/../ < 2 n - 1. 

Finally, we apply the selective phase shift to cancel the phases of A\ and A*. Defining 
<p as e %< ^ = Ai/\Ai\, we shift the phases of m basis vectors with coefficients A\ by (—<p) 
and shift the phases of m basis vectors with coefficients A\ by 0. We obtain 

RgDR e \^) = \Aor^ • • • , A 2(l+m) , ■ ■ ■ , A 2 n_ x \. (9) 
21 2m 

We write the second phase shift operator as Rg, because the phase shift angle <p depends 
on 9 and {a^}. 

If we can choose 9 to let |^4i| be equal to A , we succeed in making two sets of basis 
vectors characterized by a and a\ be weighted equally. From now, we call this series of 
operations an (RDR) operation. If we can carry out the (RDR) operations, with suitable 
parameters 9s, \n/2\ times on \ip n ), we get the uniform superposition. 

However, there are two difficulties. We can't always find a suitable 9 that lets \A±\ be 
equal to A for the (RDR) operation on an arbitrary given \ip n ). We consider the next 
lemma that shows a sufficient condition for finding a suitable 9. It gives us a hint which 
couple of sets of basis vectors do we let be weighted equally. 

Lemma 1: We define an n-qubit state \^) as 

1^) = [«o, • • • ,ai, ■ ■ ■ ,a 2 (z+m), ' ' ' > a 2"-i] for n > 2, (10) 
21 2m 

where < aj for j — 0, 1, 2(1 + m), ■ • • , 2 n — 1 and a < a\. The basis vectors of Eq. fllOD 
are G {0, 1}"}. We assume that the number of elements is equal to 21 and the 

number of elements a\ is equal to 2m, where I > 1, m > 1 and I + m < 2 n ~ 1 . We write a 
sum of all coefficients by 

2 n -l 

S = 2la + 2ma 1 + ^ aj. (11) 

j=2(l+m) 

If the following condition is satisfied, 

S - 2 n ~ 2 (a + ai ) > 0, (12) 

we can always make 2(1 + m) basis vectors whose coefficients are a® or a\ be weighted 
equally by the (RDR) operation in which R and R are applied on 2m basis vectors with 
a x . 

Proof: RgDR g \^) is given by Eq. (|]) and Eq. @. To evaluate a difference between A% 
and |^4i| 2 , we define 

f(9) = 2- 2 (^-|A 1 | 2 ) 

= (2/ao + 2ma 1 cos# + C)(aicos#-a ) - 2 n ~ 2 (a 2 - a 2 ). (13) 

If f(9) = 0, A 2 is equal to \A Y \ 2 . We estimate f(0) and /(tt/2), 

/ f(0) = (a 1 -a )[S-2 n ~ 2 (a + a 1 )}, 

1 /(tt/2) = -a (2la + C)-2 n - 2 (al-al) < 0. [ ^ } 
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If S - 2 n - 2 (a + ai) > 0, there is < 6 < (tt/2), which satisfies A 2 = \Ax\ 2 . If signs of A 
and | Ax | are different from each other, the phase shift by 7r on basis vectors with negative 
coefficients is done. ■ 
To find the suitable sequence of sets of basis vectors that we make be weighted equally, 
we take the following procedure. (For n = 2, 3, the condition of Eq. ( |T2"D is always satisfied. 
Therefore, we consider the case of n > 4.) We describe a given state \ip n ) by Eq. (H), 
where n > 4 and > for < k < \n/2\. Let a min be the minimum coefficient among 
{afc} and a min+ i be the coefficient next to a min (0 < a min < a min+ i < aj, where aj is any 
coefficient of \i/; n ) except a min and a min+ i). Because the number of different coefficients 
in {a^.} is equal to ([n/2\ + 1), it takes 0{n) steps to find a m i n and a m i n+ i on classical 
computation. 

1. If S < 2 n ~ 2 (a min + a min+ i), we get S < 2 n ~ 2 (cij + aj) for Vi, j. In this case, it can't 
be guaranteed to find a good 9 for the (RDR) operation. We take another technique 
explained in the next section. 

2. If S > 2 n ~ 2 (a min + a min+ i), we can find a good 6 for the (RDR) operation and 
get a relation, A 2 min = \A min+ i\ 2 . Because Eq. ( |T3"D is an equation of the second 
degree for cos#, we can obtain 9 with some calculations. In this case, though 
we may make other couples of sets of basis vectors be weighted equally, we neglect 
them. Shifting phases of basis vectors which have negative coefficients by 7r after the 
(RDR) operation, we obtain a state whose all coefficients are nonnegative. There 
are \n/2\ kinds of new coefficients in the state after these operations, and we can 
derive them from Eq. (H) and Eq. @ with poly(n) steps by classical computation 
(poly(n) means polynomial in n). We can check whether the condition of Lemma 1 
is satisfied or not again. 

4 The case where the sufficient condition isn't satis- 
fied 

In this section, we consider how to make a couple of sets of basis vectors be weighted 
equally in the case where the state doesn't satisfy the condition of Lemma 1, S > 
2 n ~ 2 (a min + a min+ i). We develop a technique that adjusts amplitudes of basis vectors 
and transforms the state to a state that satisfies the sufficient condition. This is an 
application of Grover's amplitude amplification process. 

For example, we consider the state which has two kinds of coefficients, 

= [ao, • • • , ax, • ■ j for n > 4, where < a < ax- (15) 
(2 n - 1) t 

We assume that the number of elements ax is equal to t, where 2 < t < 2 n — 2, and t is 
even. If < t < 2 n ~ 2 and [(3 • 2 n ~ 2 - t)/{2 n ~ 2 - t)]a < ax (ax is bigger enough than a ), 
we obtain S < 2 n ' 2 (a Q + ax) for 

In this case, applying D (we use the property of the inversion about average operation), 
and then, shifting the phase by 7r on basis vectors which have negative coefficients, we 
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Figure 1: A variation of coefficients under the (R^D) transformation. 



can reduce a difference between new coefficients, Bo and B\, as shown in Figure [I]. We 
write this phase shift operation by R^. It can be expected that [S — 2 n ~ 2 (a min + a m i n+ i)} 
gets bigger by applying (R n D) successively. The next lemma shows it clearly. 



, 02"-l 



for n > 4, 



(16) 



2 n — 1 . We assume that the number of 



Lemma 2: We consider a state, 

1^) = [ ap, - • ,ai,- ■ , a 2 (i+ m ), 
21 2m 

where < < a± < aj for j = 2(1 + m), • ■ 
elements ao is equal to 21 and the number of elements a± is equal to 2m, where / > I, 
m > 1, and I + m < 2 n ~ 1 . We also assume S, a sum of all coefficients of \^), has the 
relation, 

S -2 n ~ 2 (ao + ai ) < 0. (17) 

Applying the inversion about average operation D on\ty), and then, doing the phase shift 
transformation by 7r on basis vectors which have negative coefficients, we obtain 



iUW = [B ,---,B 1 , 



B 



2(l+m)i 



We define S as a sum of all coefficients of R^D^). We also define 



=(o) 
,(i) 



(2/ - 2 n - l )a + (2 n - 2l)a u 
(21 - 2 n - 1 )B + (2 n - 2l)B x . 



(19) 
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1. We get < B < B x < Bj for j = 2(1 + m), ■ ■ ■ , 2 n - 1 and 

[S - 2 n -\B + Bx)} -[S- 2 n ~ 2 (a + ax)] > e (0) > 0. (20) 

2. We obtain the relation, 

e« - c<°) > ^V~ 2 (a + ai ) - 5] > 0. (21) 

Proof: We can derive D\^f) = [a' , ■ ■ ■ , a[, ■ ■ ■ , a2(/+m)) ' " j a 2"-x] ? where 



(22) 



It is clear that 5 — 2™ 1 a > 0. Using the assumption of Eq. fllTD, we obtain 5 — 2™ x a fc < 
for V/c 7^ 0. Therefore, we get R n D\ty) of Eq. (ITS) , where 





= S- 


2"- 1 a , 




= S- 






= S- 


2™- 1 a i , 



2"- 1 J B = S-2 n ~ 1 a Q , 

2 n ~ 1 B 1 = -S + 2 n - 1 a 1 , (23) 
2 n - 1 B J = -S + 2 n - 1 a j , for 2(/ + m) < j < 2™ - 1. 



We can derive a difference of B\ and So with the assumption of Eq. (Jri 

2«-i( jBl - B ) = -2[S - 2 n ~ 2 (a + a x )] > 0. (24) 

It is clear that B\ < Bj for j = 2(1 + m), ■ ■ ■ , 2 n — 1. We obtain the relation, < So < 
Si < S, for j = 2(1 + m), • ■ ■ , 2" — 1. 
Since 

5 = ^r(5 -2 n " 1 a ) -S, and S + S x = Ql - a 0) 

we can derive A that is a variation of [S — 2 n_2 (ao + ax)] caused by the (R n D) operation, 

A = [S-2 n - 2 (B + B 1 )]-[S-2 n ~ 2 (a + a 1 )} 
2/ 

= 2(— - 1)5 - (4/ - 2- 1 K (25) 

To estimate A precisely, we prepare some useful relations. From the definition of S, 

we get 

2™-l 

S = 2la + 2ma x + a s > 2la + (2 n - 2l)a x . (26) 

j=2(l+m) 

Using the assumption of Eq. fllTD and Eq. (p6|) , we can derive the relation, 

> 5-2 n ~ 2 (a + ai) 

> 2/a + (T - 2l)a x - 2 n ~ 2 (a + a x ) 

= (21 - 2"~ 2 )a + (3 • 2 n ~ 2 - 2l)a x . (27) 
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We modify the relation of Eq. (p7[) and get a rougher relation, 

> 2la - 2 n ~V + (3 ■ 2 n - 2 - 2l)a x = 2la + (2 n ~ l - 2l)a x . (28) 

Because < a < cti, we obtain 21 — 2 n ~ 1 > 0. Seeing this relation and Eq. (|27| ) again, 
we also obtain 

21 > 3 ■ 2 n ~ 2 . (29) 
Here, we can estimate A. Because of Eq. (p9"D, we can substitute Eq. ( p6|) for Eq. fl25|), 

A > 2(^ T -l)[2/a + (2"-2/)a 1 ]-(4/-2"- 1 )a 

= 2^i(4Z - 3 • 2™- 1 )[(2/ - 2"" 2 )a + (3 • 2 n ~ 2 - 2/)oi] 

+2 n " 2 ( ai -a ). (30) 

Seeing Eq. (E§p, we find 3 • 2 n ~ 2 < 21 < 2 n . Therefore, we can derive the relation, 
< (4/ - 3 ■ 2 n ~ l ) < 2 n ~ l . From Eq. Q27|) and Eq. (0), we can estimate A, 

A > [(2/ - 2 n - 2 )a + (3 ■ 2"- 2 - 2l)a x ) + 2 n ~ 2 ( ai - a ) = e (0) > 0. (31) 

The first result is derived. 

From the definition of Eq. (|l^) and Eq. (^3|), (|26|) , (p7|), we can estimate the difference 
between e^ -* and 

e«-e< 0) = ^- T [(4/-3-2"- 1 ) 1 S-2"a (2/-2"- 1 )] 

> 2^i{(4Z - 3 • 2"- 1 )[2/a + (2™ - 2Z) ai ] - (21 - 2^ 1 )2"a } 

2 n — 2/ 

= -^H( 3 • 2 ^ 2 " 2/ )«i + (2i " 2 n ~ 2 )ao] 
2 n — 2/ 

> ^z^[2- 2 (a + a 1 )-5]>0. (32) 

The second result is derived. ■ 
Because of Lemma 2, doing the (R n D) transformations successively, we can make 
[S — 2 n_2 (ao + a±)] be nonnegative. We explain this matter as follows. We consider the 
state |\l/( )) specified with coefficients, < ao < a\ < ctj for 2(1 + m) < j < 2 n — 1, and 
assume S — 2 n ~ 2 (a + ai) < 0. We apply (R n D) on |\^ (0) ) and obtain described with 

coefficients, < < B[ 1] < Bf ] for 2(1 + m) < j < 2 n - 1. Because of Lemma 2.1, we 
obtain 

[SO-) _ T- 2 (B ( v ] + B[ 1] )] -[S- 2 n - 2 (a + a,)] > e< > > 0, (33) 

where is a sum of all coefficients of l 1 !^ 1 )). Then, we assume SO)-2 n - 2 (Bi 1) +B[ 1) ) < 0. 
After applying (R W D) on we get |tf< 2 )) specified by < B^ < b{ 2) < Bf for 

2(1 + m) < j < 2 n — 1. Because of Lemma 2.2, we get 

[S (2) _ 2 »-2( S (2) + B (2))] _ [S 0) _ 2 n - 2 (B^ + BP)] > e« > e<°> > 0. (34) 
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Consequently, if - 2 n ~ 2 (B^ ] + B{ k) ) < 0, [S^ - 2 n ~ 2 (B {k+X) + 5f )] increases 
by e^(> 0) at least, k stands for the number of the (R^D) transformations applied on 
the state and = S, = a , B^ = Because is defined by {a , a{\ and I, 
is a definite finite value and positive. Repeating the (R n D) finite times, we can certainly 
make [S^ - 2 n ~ 2 ( J B^ ) + B[ k) )} be nonnegative. 

From Eq. ( |2"2"| ) and Eq. during the (R n D) iteration, we find that the phase shift 
is applied on the same basis vectors. Therefore, the (R n D) iteration can be understood as 
the inversion of Grover's iteration. We use Grover's iteration for enhancing an amplitude 
of a certain basis. 

If [S {k) - 2 n - 2 (B { k) + B {k) )] comes to be nonnegative, we start to do the (RDR) 
operation again. Using the (RDR) operation and the (R n D) iteration, we can always 
transform \ip n ) to the uniform superposition. How many times do we need to apply 
(R^D) on a state to obtain the relation, - 2 n - 2 (B { k) + B [k) ) > 0? 

Estimating it, first, we introduce notations, 

e {k) = (2l-2 n - 1 )B^ k) + (2 n -2l)B (k \ 

jzik) = S {k) -2 n - 2 (B {k) + B {k) ). (35) 
Because of Lemma 2, if T^ k ' < 0, we get relations, 

e (*+i) _ e (*) > _[( 2 « _ 21) j2 n - 2 \T {k) > -F {k) /2 n ~ 3 > 0, (36) 

where we give the minimum of (2 n — 21) by 2. Consequently, if < J 7 ^ < ■ ■ ■ < 
< jz{K) < we es ti ma te e^ k ' (k = 1,2, ■ ■ ■) recurrently, 

e (D > _ x jM + e (f>) > _^W(>0), 

e (2) > - x ^ +e ^>-2x^>-2x^ 2 \>0), 

e (K) > - x F^-V + e ^>-Kx^ K -^>-Kx^ K \>0), (37) 

where x = l/2 n ~ 3 . From these relations, assuming Kx < 1, we obtain 

(0>)^ 2 ) > ^ + e^>(l-x)^>(l-x)^°\ 
(0>)^ > ^ 2 ~> + e^>(l-2x)^>(l-x)(l-2x)^°\ 

jriK+l) > T {K) + e (K) > (1 _ Kx)T {K) > _ kx)J*°l (38) 

fc=l 

If J z( - K + l ) -f > 0, we obtain J r ( K + 2 ) > and we can conclude we need to apply the 
(R n D) transformation (K + 2) times at most. To derive the upper bound on times we 
have to apply the (R n D) transformations, we estimate and J 7 ^ ' , 

e (0) > ( 2 n _ 2 ^ fli > 2ai) 

J^ (0) > 2 n a - 2 n - 2 (a + ai) > -2 n ~ 2 ai , (39) 
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and we obtain 

K K 

jriK+i) + e (o) > JJ(i _ kx)^ 0) + e (0) > -2 n - 2 ai [n(l - kx) - x). (40) 

k=l k=l 

Therefore, to estimate the lower bound on K for p( K+l ) + > 0, we have to derive 
the lower bound on K for the large n (small x) limit, where 

K 

Y[(l-kx)<x for0<x<l. (41) 

k=i 

Because lim m _^ +00 [l — (l/m)]~ m = e(> 2), if xq is small enough, we obtain 



^jrv-M<( 1 - T7 ^)^<i 



for < Vx < Xq <C 1 (|yl/x] is the minimum integer that does not below Jl/x). 
Remembering x = l/2 n ~ 3 , we get 

x>[ n C 1 - kx )T' 3 > II C 1 - (43) 

k=\Jl/i\ k=1 



Consequently, the lower bound on K is (n — 2) \yl/x~\ — 1 ~ 0(n2 n//2 ). We have to apply 
the (-Rvr-D) transformation 0(n2 n / 2 ) times at most. (See §Appendix B.) 

Using Eq. (|23|) , we can compute {-Bfc} with poly(n) steps by classical computation, 
because the number of different coefficients in {Bk} is equal to ([n/2\ + 1). 



5 The whole procedure 

In this section, we show the whole procedure for building \%p n ) and give a sketch of 
implementation for our procedure. We also show that we can use it for building more 
general entangled states. 

As a result of discussion we have had, we obtain the whole procedure to build \ip n ) 
as follows. (We describe the procedure reversely from \ip n ) to the uniform superposition. 
Throughout our procedure, we take G {0, 1}™} as basis vectors.) 

1. We consider an n-qubit register that is in the state of \tp n ) for an initial state. (We 
assume all coefficients of basis vectors are positive or equal to 0.) 

2. If the state of the register is equal to the uniform superposition, stop operations. If 
it is not equal to the uniform superposition, go to step 3. 

3. Let a m i n be the minimum coefficient for basis vectors in the state of the register and 
a min+i be the coefficient next to a min . Examine whether a min and a min+ i satisfy 
the sufficient condition of Lemma 1 or not. If they satisfy it, carry out the (RDR) 
operation, shift the phases of basis vectors which have negative coefficients by ir, 
and then go to step 2. If they do not satisfy it, go to step 4. 
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4. Apply the (R^D) transformation on the register and go to step 3. 

Before executing this procedure, we need to trace a variation of coefficients of basis 
in each step by classical computation, because we have to know which basis vectors have 
the coefficients a min and a m i n+ i, find the phase shift parameter of the (RDR) operation, 
and so on. From these results, we construct a network of quantum gates. The amount 
of classical computation is comparable with the number of steps for the whole quantum 
transformations. 

We now sketch out the points of networks of quantum gates for our procedure. Because 
it is a chain of phase shift transformations and Grover's operation Ds, we discuss the 
networks of quantum gates for them. 

First, we discuss the phase shift transformation. In the (RDR) operation, we shift the 
phases by 9 on half of basis vectors which have coefficients a k (as a m i n+ i) and by (—9) 
on the other half of them. Constructing networks for Rg, we prepare two registers and a 
unitary transformation Uf, 

\x) ® \y) \x) ® \y © f(x)), (44) 

where the first (main) register is made from n qubits, the second (auxiliary) register is 
made from m = |Tog 2 (n + 1)] qubits initialized to |0 • • • 0), and 

f(x) = (the number of "1" in the binary string of x). (45) 

Obtaining f(x) on classical computation, we need 0(nm) ~ 0(nlog 2 n) classical gates 
(XOR, and so on) and 0(m) ~ 0{\og 2 n) other auxiliary classical bits. Therefore, 
we can construct Uf with 0(n\og 2 n) elementary quantum gates ( |] (14| |15| and see 
§ Appendix A). 

To execute the selective phase shift efficiently, we apply it on the second register 
instead of the first register. Because the phase shift matrix defined in Eq. @ is diagonal, 
we can do this way. After shifting the phases, we apply U / again and initialize the second 
register. Unnecessary entanglement between the first and the second register is removed. 

To see these operations precisely, we apply Uf on \ip n ) defined in Eq. (0). We get 

\n/2\ 

Wn)®|o) = u f Y. a h \k) s ® |0) 

k=0 

EiT 1)/2 a k [\k) ®\k) + \n-k)®\n-k}), {n is odd) 
Y,tS~ l aj\k) ®\k) + \n-k)®\n- k)) 

+a n / 2 \n/2) ® |n/2), (n is even) 

where \k) is an equally weighted superposition of k excited qubits {\k) s = \k) + \n — k) 
except for \n/2) s = \n/2) where n is even). We shift the phases of basis vectors \k), 
\n — k) on the second register, instead of \k), \n — k) (that contain 2(£) binary states) 
on the first register (where k ^ n/2). This implementation reduces the number of basis 
vectors on which we apply the phase shift operation from 2(^) to 2, and we can save 
elementary quantum gates. Using another auxiliary qubit, we can carry out the phase 
shift with 0(log 2 n) elementary quantum gates ([0[|16] and see §Appendix A). If n is 
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even and k = n/2, we can't decide which basis vectors we have to shift the phases by 9 
or (—6). In this case, we refer to not only \n/2) on the second register but also the first 
qubit of the first register (cf. I^) defined in Eq. ([]])). 

Then, we discuss how to construct the quantum network of D. It is known that D 
can be decomposed to the form, D = —WRW, where W = (g) • • • ® H^ n > (the Walsh- 
Hadamard transformation on n qubits of the main register), and R is a phase shift by it 
on 1 • • - 0) of n qubits||. D takes 0{n) steps. 

We repeat the (R^D) transformation 0(n2 n / 2 ) times at most before the (RDR) oper- 
ation. If we do the (R^D) iteration before every (RDR), we carry out it [n/2\ times. 
Therefore, the (R^D) iterations take the main part of the whole steps. Because (R n D) 
takes 0(nlog 2 n) steps, we need 0((n 3 log 2 n) x 2 n//2 ) steps for the whole procedure in 
total at most. 

Finally, we consider the case where all of {a*,} are neither positive nor real. Doing 
the selective phase shift on the basis vectors with complex or negative real coefficients in 
\ip n ) to cancel the phases, we obtain a superposition whose all coefficients are real and 
nonnegative. After this operation, we can apply our procedure on the state. 

In our method, we don't fully use the symmetry of \ip n )- Essential points that we use 
are as follows. First, the number of basis vectors that have the same coefficient is always 
even. Second, we can efficiently shift the phase of half the basis vectors that have the same 
coefficients. Third, the number of different coefficients {a^} is poly(n). Therefore, we can 
apply our method to build more general entangled states that have above properties. 



X={0,l} n f Y={0,l} m 




Figure 2: A function / defined in Eq. (|4"T|). 



Here, we discuss applying our method for building more general entangled states than 
\ip n )- We consider an entangled state defined by a function /, as follows, 

/: A = {0,1}" = {0,l} m , (46) 

where m = |~log 2 (M+ 1)] + 1 and M is polynomial in n. We assume we can label elements 
of image caused by / from A = {0, 1}™ by {(0, ±), (1, ±), • ■ • , (M, ±)}. We also assume 
the number of A's elements mapped to (k, +) and the number of them mapped to (A;, — ) 
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are equal to Ik for k — 0, • • ■ , M, where 2 J2k=o = 2 n . We can describe the function / 
by 

f(x(k,e,()) = (k,e), (47) 

where k = 0, 1, • • • , M, and e = ±, and £ = 1, • • • , as shown in Figure ^|. 
Then we consider the following n-qubit partially entangled state, 

M h 

W-EEE c #- e .0> forn>2, (48) 

fc=0 e=± C=l 

where {ck} are complex. The number of sets of basis vectors classified by {ck} is (M+ 1), 
and the number of basis vectors that have the coefficient Ck is 21^. 

Executing the selective phase shift efficiently, we apply Uf of Eq. (fl?| ) to write (fc, e) 
on the m-qubit second register and apply phase shift transformation on it. We can shift 
the phase by 9 or (—9) according to e in the (RDR) operation. To transform |\I/ n ) to the 
uniform superposition, we have to do the (RDR) operation M times. Consequently, the 
(R^D) transformation is repeated M x 0(n2 n l 2 ) times at most. It is desirable that M is 
poly(n). 



6 Discussion 

It is known that any unitary transformation U (G U(2 n )) can be constructed from 
0(n 3 2 2n ) elementary gates at most|14|. In comparison with this most general case, our 



method is efficient, although the number of gates increases exponentially in n. 

C. H. Bennett et al. discuss transmitting classical information via quantum noisy 
channels|8|]. It is shown when two transmissions of the two-Pauli channel are used, the 
optimal states for transmitting classical information are partially entangled states of two 
qubits. Therefore, we can expect our method is available for quantum communication. 

Grover's algorithm was proposed as a solution of the SAT (satisfiability) problems. It 
finds a certain combination from all of 2 n possible combinations of n binary variables. 
From a different view, what Grover's method does is enhancing an amplitude of a certain 
basis vector specified with an oracle for a superposition of 2 n basis vectors. In our method, 
we use Grover's method for adjusting amplitudes of basis vectors. 

We can't show whether our procedure is optimal or not in view of the number of 
elementary gates. Because we don't use the symmetry of \ip n ) enough, our method seems 
to waste steps. 

Recently, constructing approximately an optimal state for Ramsey spectroscopy by 
spin squeezing has been proposed flr?fl. This state also has symmetry like Eq. (0), and it 
is characterized by one parameter. 
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A Networks of quantum gates 

We construct networks of quantum gates for our method concretely. For notations of 
networks and quantum gates, we refer to A. Barenco et al.[H. 



A.l The network of U 



f 



Uf (defined in Eq. (|-f), (fig) or Eq. ([47])) is given by the controlled gate, which causes 
the unitary transformation on the second register under the value of the first register. 
Constructing the controlled gate of Uf with poly(n) quantum elementary gates, we can 
use our method efficiently 

We consider a network for Uf defined in Eq. ( fP]) and Eq. (f45|). f(x) represents the 
number of "1" bit in the binary string x. Writing the first (main) and second (auxiliary) 
register by \X n , X n _i, • • • , X 2 , X\) <%>\S), where \S) is made up of m = |~log 2 (n + 1)] qubits 
and initialized to [0 • • • 0), we can write the quantum networks as the following program. 
(For the notation of the program, we referred to Cleve et al.[ 

Program adder- 1 
quantum registers: 

Xi, X 2 , ■ ■ ■ , X n : qubit registers 
S: an m-qubit register 
for k = 1 to n do 

S <- (S + X k ) mod 2 m . 



rTog2(n+iyi 

qubits 
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Figure 3: The network of the adder-2 for m = [log 2 (n + 1)] = 4. 



To write a program for the addition of X k in the adder- 1, we describe the qubits of 
the second register by [SVrt-i, • • • ,Si, So), introduce other auxiliary qubits |C m _i, • • • , Ci), 
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and use Cj as a carry bit of addition at the (j — l)th bit. We can write the program as 
follows. 

Program adder- 2 
quantum registers: 

X k : a qubit register 

So, Si, ■ ■ ■ , S m -i. qubit registers 

Ci, C2, ■ ■ ■ , C TO _i: auxiliary qubit registers (initialized and finalized to 0) 
for j ' = 2 to m — 1 do 

Cj «- © (Cj-i a a,--!) 

for j = m — 1 down to 2 do 

<— Sj © Cj 
Cj - C, © (Cml A 
Si ^ Si © Ci 
C 1 ^C 1 ©(S , AX fe ) 
So ^ S © 



Because we don't use interference, we can describe these operations with a higher 
level language of classical computation. In this program, to avoid obtaining unnecessary 
entanglement, we initialize and finalize all auxiliary qubits {|Cj)} t° |0). Figure |3] shows 
a network of this program for m — 4. Repeating the quantum network of the adder-2 for 
each X k (k — 1, ■ • • , n), we can construct the adder- 1. 

We estimate the number of quantum elementary gates to construct the adder-1. In 
Figure 0, we use 2([log 2 (n+ 1)] — 1) Toffoli gates (that maps \x,y,z) \x, y, z® (x Ay))) 
and [log 2 (n +1)] controlled-NOT gates for the adder-2. Because we repeat the adder-2 
n times, the number of the whole steps for the adder-1 is equal to n{?>\\og 2 {n + 1)] — 2). 



A. 2 Construction of f\ n (R z (a)) 

From now on, we often use a /\ n (R z (a)) gate, where R z (a) is given in the form, 



R z {ol) = exp(iaa z /2) 



exp(ia/2) 

exp(— ia/2) 



(49) 



(We describe the controlled™-?/ by f\ m {U), where WU G U(2). f\ m {U) has an m-qubit 
control subsystem and a one-qubit target subsystem. It works as follows. If all m qubits 
of control subsystem are equal to |1), f\ m {U) applies U on a target qubit. Otherwise 
f\ m {U) does nothing. We can write the Toffoli gate by ^(o'x), the controlled-NOT gate 
by Ai(°"x), and any U{2) gate for one qubit by AoO Here, we consider how to construct 
it from elementary gates. 
At first, using relations, 

R z {a/2)a x R z {-a/2)a x = R z {a), and R z {a/2)R z {-a/2) = I, 

we can decompose a f\ n (R z (a)) gate to a f\ l (R z (a/2)) gate, a /\ l (R z (—a/2)) gate and two 
An-i(o"x) gates, as shown in Figure [|. Seeing Figure [|, we can decompose a Ai(-Rz(/3)) 
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n qubits [ 



R z (oc) 




Figure 4: Decomposition of a /\ n (R z (a)) gate. 



Rz(P) 



R z (p/2) R z (-p/2) I-® 



Figure 5: Decomposition of a Ai(-R 2 (/9)) gate. 



gate to an R z (j3/2) gate, an R z (—/3/2) gate and two controlled-NOT gates. We have to 
only consider how to make a f\ n -i{(J x ) gate from elementary gates on an (n + l)-qubit 
network. Especially we pay attention to the fact that there is a qubit which is not used 
by the f\ n _i(a x ) gate on the network. 

It is shown that, on an (n + l)-qubit network, where n > 6, a An-iO^) gate can be 
decomposed to 8(n — 4) Toffoli gates|T4J]. Consequently, on the (n+ l)-qubit network(n > 
6), a /\ n (R z (a)) gate can be decomposed to 16(n — 4) Toffoli gates, four f\i(a x ) gates and 
four Ao gates. Therefore, f\ n (R z (a)) takes 8(2n — 7) quantum elementary gates in total. 



A. 3 The phase shift on certain basis vectors 

Figure § shows a quantum network for the selective phase shift by 9 on a certain basis 
vector of the second register defined in Eq. fl4~4]). In Figure [| we use a f\ m (R z {26)) gate. 
Setting the auxiliary qubit being |0), /\ m (R z (29)) generates an eigenvalue exp(z#) if and 
only if the second register is in the state |1 ■ • ■ 1). This technique is called "kick back' |16| . 




Figure 6: The network of the phase shift on the second register. 



In Figure ^ a shaded box stands for the NOT-gate given by a x (|0) — > |1), |1) — > |0)) 
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or the identity transformation. Deciding which gates are set in each shaded box, a x or J, 
we can select a basis vector on which we shift the phase. 

In case m > 6, it has been already shown that a /\ m (R z (29)) gate can be constructed 
from 8 (2m — 7) quantum elementary gates at most. Seeing Figure |6|, we find that the 
selective phase shift on the second register takes 2m + 8(2m — 7) = 2(9m — 28) ~ 0(m) 
gates at most. Building \ip n ), we can carry out the phase shift on certain basis vectors on 
the second register with 0(log 2 n) steps. 

A. 4 The network of D 

Figure |7] shows a network of D. Since this network consists of 4n Ao gates and a 
A n (-R z (27r)) gate, it takes 4(5n — 14) elementary gates, in case n > 6. Therefore, D 
takes 0{n) steps. 
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Figure 7: The network of D. 



A. 5 Estimation of steps 

How many elementary gates do we need to construct \ip n ) defined in Eq. (Q) or |\l/ n ) 
defined in Eq. ( f48"D from the uniform superposition? If M (the number of sets of basis 
vectors classified by their coefficients) is poly(n), and if the function Uf defined in Eq. ( |44| ) 
can be constructed from poly(n) elementary gates, the (R n D) iterations take the main 
part of the whole steps. 

In the (R n D) transformation, we do the following operations. Applying D on the 
n-qubit first register, preparing the initialized m-qubit second register, we apply Uf on 
both of the registers as Eq. (f|4|). Then, we shift the phases of basis vectors on the second 
register. Finally, we apply Uf again to initialize the second register. It has been already 
shown D takes 4(5n — 14) steps, where n > 6. The number of steps that a network of Uf 
takes depends on the function /. For instance, when we build \ip n ), Uf needs 0(nlog 2 n) 
steps. 

Figure |8| shows a network of the phase shift by ir on the second register with negative 
coefficients for R n , in the case of building \ip n ). To inverse signs of negative coefficients, 
we shift the phase at most for [n/2\ sets of basis vectors characterized by coefficients. 
Therefore, we shift the phases of 2|_n/2j basis vectors of the second register at most. As 
a result, the network can be constructed from (2|_n/2j + 1) • |~log 2 (n + 1)] Ao gates and 
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Figure 8: The network of the phase shift by 7r on the second register for R n 



2\n/2\ Ariog (n+i)l(-^z(27r)) gates. We can carry out R n with 0(nlog 2 n) steps. Similarly, 
in the case of |^ n ), we can estimate R^ takes 0(Mlog 2 M) steps. 

Building \tp n ), we repeat the (R^D) transformation 0(n2 n / 2 )times at most before the 
(RDR) operation. If we do the (R n D) iteration before every (RDR), we carry out it 
|_n/2j times. Consequently, we need 0((n 3 log 2 n) x 2 n//2 ) steps for the whole procedure 
in total at most. 

B A variation of coefficients during the (R n D) itera- 
tion 

We explicitly derive a variation of coefficients during the (R n D) iteration for the case 
described in (|T5|), and estimate how many times do we need to apply (R n D) to make 
[SW - 2 n ~ 2 (B$ + B[ k) )} be nonnegative. We find it takes 0{2 n > 2 ) times, in spite of the 
results, 0{n2 n l 2 ) times, in §f|. 

Applying (R n D) on |^) defined by (pj), we obtain R n D\^) 
where 

2 n ~ 1 B^ ) = S - 2 n - 1 a = (2™- 1 - t)a + ta u 
2 n-i jB (i) = _ S + 2 ^ ai = -{2 n - t)a + (2 n " ] 

Referring to 0, we write t as 

t 2 n — t 
sin 2 9 = —, (cos 2 9 = - -), 

On ' v n n 

where < 9 < (vr/2), and write {ao,ai} as 

sin a cos a 



t)ai. 



(50) 



(51) 



(52) 
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where < a < (tt/2). Using (|50|), (pl|), and (|52]) , we can describe {S^,^} by 

f Bj 1 ) = (l/v/2")[cos20(sin a/ cos 9) + 2 sin cos a] = sin(a + 2fl)/ v / 2 7rr t, , * 



5} ; = (l/ v / 2 ir )[-2cos#sina + cos2#(cosa/sin#)] = cos(a + 29)/^. 



Writing coefficients of the state on which (R n D) has been applied k times as and 
Bi, we obtain 

^ fc) = sin(a + 2fcfl) ( fc) = cos(a + 2ifc6>) for A; = 0, 1, 2, • • • , (54) 

where = a , b[ 0) = a x . 

Defining = (2 n - t)B { k) + tB[ k \ we can derive 

S {k) _ 2»-2( S W + 

= (3 ■ 2 n - 2 - t)5j fc) + (t - 2 n - 2 )B[ k) 

2^{sin[a + (2fc + 1)0] - - — cos[a + (2k - 1)9}} 

2 sin 29 



,/On-2 

X " -F« (55) 



sin 20 



where 

= cos[a + (2ife + 3)0]. (56) 

Since < < (vr/2) and sin 29 > 0, it depends on a sign of whether [S^ - 
is negative or not. (With some calculations, we can confirm that (|54D 
and ( |55| ) satisfy Lemma 2.) 

Because of < a < (tt/2), if (2k + 3)9 = (tt/2), it is always accomplished that < 
and — 2 n ~ 2 (B ( Q S> + B^) > 0. Therefore, the number of times we need to apply (R^D) 
doesn't exceed kuAx-, which is given as 

k MA x = ^(l-m- (57) 



On the other hand, we can write 9 as sin0 = Jt/2 n from (|51"D, and the minimum value 
of t is 2. Taking the limit that t ~ 0(1) and n is large enough, we obtain a relation, 

7T / 2™ 

A^max ~ ~ 0(2™/ 2 ). (58) 



sin0 ~ 6> ~ v/t/2 n and 



The (R n D) transformation is repeated 0(2 n ^ 2 ) times at most to make [S^ — 2 n 2 (B^ + 
B^)} be nonnegative. 
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